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1.  INTRODUCTION 


Density  functional  theory  for  stationary  states  or  ensembles  is  a  for¬ 
mulation  of  many-body  theory  in  terms  of  the  particle  density,  p(r). 
Based  on  the  work  of  Hohenberg  and  Kohn1  (HK)  and  Kohn  and 
Sham2  (KS),  it  is  now  a  mature  subject  which  has  had  many  success¬ 
ful  applications.  In  retrospect,  the  stationary  Thomas-Fermi  theory3,4 
and  the  Hartree5  and  Hartree-Fock-Slater6  equations  can  be  viewed  as 
approximations  of  density  functional  theory. 

Time-dependent  density  functional  theory  as  a  complete  formalism7 
is  of  more  recent  origin,  although  a  time-dependent  version  of  Thomas- 
Fermi  theory  was  proposed  as  early  as  1933  by  Bloch.®  The  time- 
dependent  Thomas-Fermi  model  provides  good  first  estimates  for  a 
variety  of  physical  processes  such  as  atomic9-11  and  nuclear12'13  col¬ 
lisions,  nuclear  giant  resonances,14-16  atomic  photoabsorption,17-19 
and  plasmons  in  inhomogeneous  media.20-25  The  time-dependent 
Hartree  and  Hartree-Fock  equations26  have  been  used  extensively  for 
more  quantitative  studies  of  nuclear  dynamics27,28  and  atomic29-32 
as  well  as  molecular33  collision  processes.  However  these  methods  ne¬ 
glect  correlation  effects.  Inclusion  of  such  effects  would  require  the 
time-dependent  analog  of  the  stationary  KS  equations. 

The  first,  and  rather  successful,  steps  towards  a  time-dependent  KS 
theory  were  taken  by  Peuckert34  and  by  Zangwill  and  Soven.35  These 
authors  treated  the  linear  density  response  of  r are-gas  atoms  to  a 
time-depen  dent  external  potential  as  the  response  of  non-interacting 
electrons  to  an  effective  time-dependent  potential.  In  analogy  to  sta¬ 
tionary  KS  theory,  this  effective  potential  was  assumed  to  contain  an 
exchange-correlation  part,  Vze(rt),  in  addition  to  the  time-dependent 
external  and  Hartree  terms: 

ve//(rt)  =  v(rt)+/  dV+Vxcfrt). 


Peuckert  suggested  an  iterative  scheme  for  the  calculation  of  vxe, 
while  Zangwill  and  Soven  adopted  the  functional  form  of  the  static 
exchange-correlation  potential, 


vic(rt)  = 


8Exc\p{'t>i) 

Sp[rt)' 


where  Exe[p\  is  the  exchange-correlation  energy  functional  of  ordinary 
density  functional  theory.  For  this  functional,  Zangwill  and  Soven  em¬ 
ployed  the  local  density  approximation, 


Exe[p{Tlt)\  =  f  eXc{p{T,t))p(T,t)d?rl, 

where  tZc  [p]  is  the  exchange-correlation  energy  per  particle  of  a  static 
uniform  electron  gas  of  density  p.  The  static  approximation  is  obvi¬ 
ously  valid  only  if  the  time-dependence  of  p(rt)  is  sufficiently  slow.  In 
practice,  however,  it  gave  quite  good  results  even  for  the  case  of  rather 
rapid  time-dependence. 

Significant  steps  towards  a  rigorous  foundation  of  time-dependent 
density  functional  theory  were  taken  by  Deb  and  Ghosh36-39  and  by 
Bartolotti40-43  who  formulated  and  explored  Hohenberg-Kohn  and 
Kohn-Sham  type  theorems  for  the  time-dependent  density.  Each  of 
these  derivations,  however,  was  restricted  to  a  rather  narrow  set  of 
allowable  time-dependent  potentials  (to  potentials  periodic  in  time 
in  the  theorems  of  Deb  and  Ghosh,  and  to  adiabatic  processes  in  the 
work  of  Bartolotti). 

Finally,  a  general  formulation  covering  essentially  all  time-dependent 
potentials  of  interest  was  given  by  Runge  and  Gross.7  A  novel  fea¬ 
ture  of  this  formalism,  not  present  in  ground-state  density  functional 
theory,  is  the  dependence  of  the  respective  density  functionals  on  the 
initial  (many-particle)  state  ♦(to).  A  detailed  description  of  the  time- 
dependent  density-functional  formalism  will  be  presented  in  section 
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2.  The  central  result  is  a  set  of  time-dependent  KS  equations  which 
are  structurally  similar  to  the  time-dependent  Hartree  equations  but 
include  (in  principle  exactly)  all  many-body  correlations  through  a 
local  time-dependent  exchange  correlation  potential. 

To  date,  most  applications  of  the  formalism  fall  in  the  regime  of  lin¬ 
ear  response.  In  section  3,  we  shall  describe  the  linear-response  limit 
of  time-dependent  density  functional  theory  along  with  applications 
to  the  photo-response  of  atoms,  molecules  and  metallic  surfaces. 

Beyond  the  regime  of  linear  response,  the  description  of  atomic  and 
nuclear  collision  processes  appears  to  be  a  promising  field  of  applica¬ 
tion  where  the  time-dependent  KS  scheme  could  serve  as  an  econom¬ 
ical  alternative  to  time-dependent  configuration-interaction  calcula¬ 
tions.44-49  So  far,  only  a  simplified  version  of  the  time-dependent  KS 
scheme  has  been  implemented  in  this  context.50-53 

Another  possible  application  beyond  the  regime  of  linear  response 
is  the  calculation  of  atomic  multiphoton  ionization  which,  in  the  case 
of  hydrogen,  has  recently  been  found54,55  to  exhibit  chaotic  behaviour. 
A  full-scale  numerical  solution  of  the  time-dependent  Schrodinger 
equation  for  a  hydrogen  atom  placed  in  strong  time-dependent  elec¬ 
tric  fields  has  recently  been  reported.56  A  time-dependent  Hartree- 
Fock  calculation  has  been  achieved  for  the  multiphoton  ionization 
of  helium.57  For  heavier  atoms  an  analogous  solution  of  the  time- 
dependent  Kohn-Sham  equations  offers  itself  as  a  promising  applica¬ 
tion  of  time-dependent  density  functional  theory. 

2.  FORMAL  FRAMEWORK 

2.1.  One-to-one  mapping  between  time-dependent  potentials 
and  time-dependent  densities 

Density  functional  theory  is  based  on  the  existence  of  an  exact 
mapping  between  densities  and  external  potentials.  In  the  ground- 
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state  formalism,1  the  existence  proof  relies  on  the  Rayleigh-Ritz  min¬ 
imum  principle  for  the  er  'rgy.  Straightforward  extension  to  the  time- 
dependent  domain  is  not  possible  since  a  minimum  principle  is  not 
available  in  this  case.  The  proof  given  by  Runge  and  Gross'  for  time- 
dependent  systems  is  based  directly  on  the  Schrodinger  equation* 

<  &»(<)  -  A(t)*(t).  (i) 

We  shall  investigate  the  densities  p(r t)  of  electronic  systems  evolving 
from  a  fixed  initial  (many-particle)  state 

*(*o)  =  *0  (2) 

under  the  influence  of  different  external  potentials  of  the  form 

V{t)  =  £,r=n/<*3r  tfc(rMrt)tMr).  (3) 

In  the  following  discussion,  the  initial  time  to  is  assumed  to  be  finite 
and  the  potentials  are  required  to  be  expandable  in  a  Taylor  series 
about  to: 

«(«)  =  sjio  <>>(•■)((  -  «0)J-  w 

No  further  assumptions  concerning  the  size  of  the  radius  of  conver¬ 
gence  are  made.  It  is  sufficient  that  the  radius  of  convergence  is 
greater  than  zero.  The  initial  state  ¥o  is  not  required  to  be  the 
ground  state  or  some  other  eigenstate  of  the  initial  potential  v(rto)  = 
vo(r).  This  means  that  the  case  of  sudden  switching  is  included  in 
the  formalism.  On  the  other  hand,  potentials  that  are  switched- 
on  adiabatically  from  to  =  -oo  are  excluded  by  the  condition  (4). 

*  Atomic  units  are  used  throughout. 


-6- 


Besides  an  external  potential  of  the  form  (3),  the  Hamiltonian  in 
Eq.  (1)  contains  the  kinetic  energy  of  the  electrons  and  their  mutual 
Coulomb  repulsion: 

H[t)  =  T-hU-hV(t)  (5) 

with 

i  =  E *.nfd3'  &»(-  (6) 

and 

u  =  lZ„,rfd3rfd3r’  tirtt+V)  m 

The  number  of  electrons,  N,  is  fixed. 

We  shall  demonstrate  in  the  following  that  the  densities  p(rt) 
and  p'(rt)  evolving  from  a  common  initial  state  ¥o  under  the  influ¬ 
ence  cf  two  potentials  v(rt)  and  */(rt)  are  always  different  provided 
that  the  potentials  differ  by  more  than  a  purely  time-dependent  (r- 
independent)  function:! 

v(rt)  ±  v'(rt)  +  c(i).  (8) 

Since  both  potentials  can  be  expanded  in  a  Taylor  series  and  since  they 
differ  by  more  than  a  time-dependent  function,  some  of  the  expansion 
coefficients  vy(r)  and  v'(r)  (cf.  Eq.  (4))  must  differ  by  more  than 

t  If  v  and  v'  differ  by  a  purely  time-dependent  function,  the  resulting 
wave  functions  ¥(t)  and  ♦/(t)  differ  by  a  purely  time-dependent  phase 
factor  and,  consequently,  the  resulting  densities  p  and  p1  are  identical. 
This  trivial  case  is  excluded  by  the  condition  (8),  in  analogy  to  the 
ground-state  formalism  where  the  potentials  are  required  to  differ  by 
more  than  a  constant. 
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a  constant.  Hence  there  exists  some  smallest  integer  fc>0  such  that 

vjfc(r)-ojfc(r)  =  fp£  [»(rt)-V(rt)lt«*o  #  const.  (9) 

From  this  inequality  we  shall  deduce  in  a  first  step  that  the  current 
densities  j(rt)  and  j'(rt)  corresponding  to  u(rt)  and  v'(rt)  are  differ¬ 
ent.  In  a  second  step  it  will  be  demonstrated  that  the  densities  p(rt) 
and  p'{vt)  are  different. 

Since  we  consider,  by  construction,  only  states  ¥(t)  and  ^'(t)  that 
evolve  from  the  same  initial  state  'to*  the  current  densities  j  and  j' 
as  well  as  the  densities  p  and  p1  are  identical  at  the  initial  time  t0: 


j(r*o)  =  j'(r*0)  =  <  'toli(r)!*0  >  =  Jo(r),  (10) 

P(rt0)  =  p'(rt0)  =  <  't0|p(r)i't0  >  =  p0(r).  (11) 

j  (r)  denotes  the  paramagnetic  current-density  operator 

J(r)  =  -i  E»{*5?(r)lVA,(r)|-[V*+(r)l^(r)},  (12) 

and  p(r)  is  the  usual  density  operator 

P(r)  =  U3) 

The  time  evolution  of  the  current  densities  j  and  j;  is  most  easily  com¬ 
pared  by  means  of  the  equations  of  motion 

jfJ(rt)=  -.<*(l)lli(r),H(l)||*(l)>  (14) 

JfJV)  =  -.<*'(<)|U(r).H'(()]|*'(l)>.  (15) 
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Taking  the  difference  of  these  two  equations  at  the  initial  time  to  one 
finds 

Jfli(r<)-j'(rl)lt-to  =  -*<*olli(r),(H(io)-H'(<o))]l*o> 

=  -<>o(r)V(ti(rto)-»'(«o)).  (16) 

If  the  potentials  u(r t)  and  v'(rt)  differ  at  t  =  <o  (i.e.,  if  Eq.  (9)  is  satis¬ 
fied  for  k  =  0)  then  th«*  right-hand  side  of  Eq.  (16)  cannot  vanish  iden¬ 
tically.  Consequently  ^rt)  and  j'(r t)  will  become  different  infinitesi¬ 
mally  later  than  to-  If  the  minimum  integer  k  for  which  Eq.  (9)  holds 
is  greater  than  zero,  one  has  to  evaluate  the  (fc+l)th  time-derivative  of 
j  and  j'  at  f0.  To  this  end,  the  quantum  mechanical  equation  of  motion 

i  <  «(i)|6(i)|*(i)  >  =  <  *(t)l£d(t)-<[d(«),  j(!)n»(o  > 

is  applied  ( k  +  1)  times;  the  first  time  with  Oi(t)  =  j(r),  the  second 
time  with  02(t)  =  -t[j(r), //(*)],  etc.  After  some  straightforward  al¬ 
gebra  one  obtains 

[j(rf)-j'(rt)]<=tc 

=  ¥■  0.  (17) 

where  the  last  inequality  follows  from  (9)  Again  the  inequality  (17) 
implies  that  j(rt)  and  j'(rt)  will  become  different  infinitesimally  later 
than  to-  This  completes  the  proof  for  the  current  densities. 

In  order  to  prove  that  the  densities  p(rt)  and  p'(r t)  are  different, 
one  employs  the  continuity  equation: 


$(,(«)-✓(«)!  =  -VP(«)-j'(rt)|. 


(18) 
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Evaluating  the  (fc-f 1)  th  time-derivative  of  Eq.  (18)  at  the  initial  time 
*0i  one  finds  by  insertion  of  Eq.  (17) 

[/>(«) -p'(rt)]t=to  =  V-fpoMVtMr)]  (19) 

where 

^fc(r)  =  ^f(v(r<)-v'(r<)]t=t0  =  wfc(r)~v*(r).  (20) 

It  remains  to  be  shown  that  the  right-hand  side  of  Eq.  (19)  cannot 
vanish  identically  if  the  condition  (9),  ty^(r)  ^  const,  is  satisfied.  The 
proof  is  by  reductio  ad  abaurdum:  Assume  that  the  right-hand  side  of 
Eq.  (19)  vanishes  and  consider  the  integral 

fd3r  p0(r)[Vu;fc(r)]2  (21) 

which,  by  means  of  Green’s  theorem,  can  be  written  as 

=  -fd3r  wk(r)V-lp0(r)Vwk(r)J  +  /dS-(po(r)u'*(r)Vu;it(r)).(22) 

As  long  as  one  deals  with  realistic,  i.e.  experimentally  realizable  po¬ 
tentials  v(rt),  the  surface  integral  in  (22)  vanishes:  Experimentally 
realizable  potentials  are  always  produced  by  real  (this  means,  in  par¬ 
ticular,  normalizable)  external  charge  densities  pext( rt)  so  that,  in  the 
Coulomb  gauge, 

v(rt)  =  (23) 

A  Taylor  expansion  of  this  equation  about  t  =  t o  demonstrates  that 
the  expansion  coefficients  v*(r),  and  hence  the  functions  itffc(r),  fall 
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off  at  least  as  1  /r  asymptotically  so  that  the  surface  integral  vanishes. 
(For  more  general  potentials  see  remark  (iii)  below.)  The  first  term 
of  (22),  on  the  other  hand,  is  zero  by  assumption.  Thus  the  complete 
integral  (21)  varnishes  and,  since  the  integrand  is  nonnegative,  one  con¬ 
cludes 


/>o(r)(Vt^Jfc(r)]2  =  0 

in  contradiction  to  u>*(r)  ^  const.*  This  completes  the  proof  of  the 
theorem.  A  few  remarks  are  in  order  at  this  point: 

i)  By  virtue  of  the  1-1  correspondence  established  above  (for  a 
given  *o).  the  time-dependent  density  determines  the  external  po¬ 
tential  uniquely  up  to  within  an  additive  purely  time-dependent  func¬ 
tion.  The  potential,  on  the  other  hand,  uniquely  determines  the  time- 
dependent  wave  function,  which  can  therefore  be  considered  as  a  func¬ 
tional  of  the  time-dependent  density, 


»(!)  = 


(24) 


where  ^[p]  is  unique  up  to  within  a  purely  time-dependent  phase  fac¬ 
tor.  As  a  consequence,  the  expectation  value  of  any  quantum  mechani¬ 
cal  operator  Q(t)  is  a  unique  functional  of  the  density, 


Q{t)  =  <  *M(t)l3(t)l*M(0  >  =  QMW;  (25) 


the  ambiguity  in  the  phase  of  ¥{p]  cancels  out.  As  a  particular  ex¬ 
ample,  the  right-hand  side  of  Eq.  (14)  can  be  considered  as  a  density 

*  It  is  assumed  here  that  po(r)  does  not  vanish  on  a  set  of  positive 
measure.  (Otherwise  no  contradiction  is  reached  since  po(r)  could  be 
zero  in  exactly  those  regions  of  space  where  u/fc(r)  /  const.)  Densi¬ 
ties  that  vanish  on  a  set  of  positive  measure  correspond  to  potentials 
having  infinite  barriers.  Such  potentials  have  to  be  excluded  from  the 
ordinary  ground-state  theory  as  well  (see,  e.g.,  Ref.  58). 
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functional  which  depends  parametrically  on  r  and  t : 

PM(rt)  =  -.<*M(0l[iM,«(0ll*M(<)>-  (26) 

This  implies  that  the  time-dependent  particle  and  current  densities 
can  always  be  calculated  (in  principle  exactly)  from  a  set  of  “hydrody- 
namical”  equations: 

^p(rt)  =  -V-j(rt) 

(27) 

=  P[p](rt). 

In  practice,  the  functional  P[p]  is  of  course  only  approximately  known. 

ii)  In  addition  to  their  dependence  on  the  density,  the  functionals 
Q[p],  defined  by  Eq.  (25),  implicitly  depend  on  the  initial  state  ^0- 
No  such  state  dependence  exists  in  the  ground  state  formalism.  Of 
course  one  would  prefer  to  have  functionals  of  the  density  alone  rather 
than  functionals  of  p(rt)  and  ¥<}.  It  should  be  noted,  however,  that  for 
a  large  class  of  systems,  namely  those  where  ¥o  is  a  non-degenerate 
gound  state,  Q[p\  is  indeed  a  functional  of  the  density  alone.  This  is 
because  any  non-degenerate  ground  state  ¥<)  is  a  unique  functional  of 
its  density  po(r)  by  virtue  of  the  traditional  HK  theorem.1  One  has  to 
emphasize  that  this  class  of  systems  also  contains  all  cases  of  sudden 
switching  where  ¥<>  1®  *  ground  state,  but  not  the  ground  state  of  the 
initial  potential  v(rto).f 

iii)  It  is  essential  for  the  proof  of  the  1-1  correspondence  that  the 
surface  integral  in  (22)  vanishes.  As  long  as  one  considers  only  realis¬ 
tic  physical  potentials  of  the  form  (23),  a  vanishing  surface  integral  is 

t  Whether  am  exclusive  dependence  on  the  density  alone  can  be  demon¬ 
strated  for  an  even  larger  class  of  systems  is  am  open  question  at 
present. 
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guaranteed.  However,  if  one  allows  more  general  potentials,  the  sur¬ 
face  integral  does  not  necessarily  vanish.  Consider  a  given  initial  state 
#0  that  leads  to  a  certain  asymptotic  form  of  Pq{t).  Then  one  can 
always  find  potentials  which  increase  sufficiently  steeply  in  the  asymp¬ 
totic  region  such  that  the  surface  integral  does  not  vanish.  For  such 
cases,  the  right-hand  side  of  Eq.  (19)  can  be  zero,  as  was  demonstrated 
by  Xu  and  Rajagopal59  (see  also  Dhara  and  Gosh®0)  with  the  explicit 
examples 

p0(r)  =  Ae~Xr /(Ar)2,  wk(r)  =  BeXr,  i4,A>0  (28) 


and 


Po(r)  =  const,  ttffc(r)  satisfying  V2u7j.(r)  =  0.  (29) 

If  one  intends  to  include  potentials  other  than  those  given  by  Eq.  (23), 
an  additional  condition  must  be  imposed  to  ensure  that  the  surface 
integral  in  (22)  vanishes.  For  a  giver;  .nitial  state  'to.  the  allowable  po¬ 
tentials  must  satisfy  the  condition* 

{r2Po(r)wy(r)Vvy(r)}  =  0  (30) 

for  all  coefficients  vy(r)  of  the  Taylor  expansion  (4).  Thus,  the  density 
functionals  Q[p]  given  by  Eq.  (25)  not  only  depend  explicitly  on  the 
initial  state  '9q\  the  latter  also  determines  the  set  of  allowable  poten¬ 
tials  and  hence  the  domain  of  densities  for  which  the  functionals  Q[p] 
are  defined.  One  has  to  emphasize,  however,  that  the  cases  excluded 
by  the  additional  condition  (30)  are  largely  unphysical.  The  examples 

*  The  condition  pQWkVivic-+0,  given  by  Dhara  and  Ghosh®0  for  r— *•  oo, 
is  n£t  sufficient  to  ensure  a  vanishing  surface  integral  since  it  does  not 
account  for  the  fact  that  the  spherical  surface  increment  dS  increases 
proportional  to  r2. 
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(28)  and  (29)  of  Xu  and  Rajagopal,  for  instance,  lead  to  an  infinite! 
potential  energy  / p(rt)v(r*)d3r  per  particle  in  the  vicinity  of  t  =  t0. 
It  would  be  desirable  to  prove  the  1-1  correspondence  under  a  physical 
condition,  such  as  the  requirement  of  a  finite  potential  energy  per  par¬ 
ticle,  rather  than  under  the  more  restrictive  mathematical  condition 
(30).  Whether  such  a  proof  can  be  constructed  is  currently  an  open 
question. 

2.2.  Variational  principle 

The  solution  of  the  time-dependent  Schrodinger  equation  (1)  with 
initial  condition  (2)  corresponds  to  a  stationary  point  of  the  quantum 
mechanical  action  integral 

A  =  f‘,‘  it  <  *(0I<&  -  tf(!)l*(!)  >  .  (31) 

Since  there  is  a  1-1  mapping  between  time-dependent  wave  functions, 
'l'(t),  and  time-dependent  densities,  p(rt),  the  corresponding  density 
functional 

AM  =  /£<#<  *[ol(<)l'Jf  -  £(<)I*W(I)  >  (32) 

must  have  a  stationary  point  at  the  correct  time-dependent  density 
(corresponding  to  the  Hamiltonian  H{t)  and  the  initial  state  ¥o)- 
Thus  the  correct  density  can  be  obtained  by  solving  the  Euler  equation 


with  appropriate  boundary  conditions.  The  functional  A[p ]  can  be 


f  Note  that  a  function  satisfying  V2wjfe( r)  =  0  for  all  r  cannot  vanish 
asymptotically  (except  in  the  trivial  case  w* (rj  =  0),  nor  can  it  ap¬ 
proach  a  constant  asymptotic  value  (except  in  tne  case  u/^(r)  =  const). 
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written  as 

A[p\  =  B[p\  -  dtfd3r  p(rt)v(rt)  (34) 

with  a  universal  (^-dependent)  functional  B[p],  formally  defined  as 
B\p\  =  /!‘  dt  <  »(«](<)l>£-#,-0|*M(»)  >  .  (35) 


2.3.  Time-dependent  Kohn-Sham  scheme 

A  particularly  important  application  of  the  stationary  action  prin¬ 
ciple  described  in  the  last  section  is  the  derivation  of  a  time-dependent 
KS  scheme.  To  this  end  we  first  consider  a  system  of  non-interacting 
fermions  subject  to  an  external  potential  v(rt).  The  exact  many- 
particle  wave  function  is  then  a  time-dependent  Slater  determinant 
$(0  which,  by  virtue  of  the  Runge-Gross  theorem,  is  a  functional, 
$(p](t),  °f  the  time-dependent  density.  (Note  that  the  argument  of 
section  2.1  holds  true  for  any  particle-particle  interaction  U ,  in  par¬ 
ticular  also  for  U  =  0.  The  functionals  P[p],  B[p],  etc.,  however  de¬ 
pend  on  the  particular  U  considered.)  The  action  functional  for  non¬ 
interacting  particles  then  has  the  form 

Aa[p\  =  B,[p]  -  /5‘  dtfd3r  p(rt)v(rt)  (36) 

where 

B.W  =  f‘,‘dt<  >  •  (37) 

According  to  the  variational  principle  of  section  2.2,  the  exact  density 
can  be  obtained  by  solving  the  Euler  equation 
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(38) 


Next  we  assume!  that  a  single-particle  potential  v,(rt)  exists  such  that 
the  interacting  density  is  identical  with  the  density 


P(rt)  =  £n=  1  l^n(rt)i2 


(39) 


of  the  non-interacting  system  described  by  the  Schrodinger  equation 

<Pn{*t)  =  (-^-  +  t;a(rt))v3n(rt).  (40) 

If  such  a  potential  exists  it  must  be  unique  (by  virtue  of  the  Runge- 
Gross  theorem)  and,  .  view  of  Eq.  (38),  it  can  be  represented  as 

Mr*)  =  I  £(rt)=p(rt),  (41) 

where  the  functional  derivative  of  Ba  has  to  be  evaluated  at  the  exact 
interacting  (=  non-interacting)  density  p(rt). 

To  determine  the  potential  vs  (rt)  more  explicitly  one  employs  the 
variational  principle  for  the  interacting  system.  First  one  rewrites  the 
total  action  functional  (34)  as 

A[p]  =  B,[p\  -  /jj  dtfd3rp(rt)v(rt) 

-ifi ;  elskgja (42) 

where  the  “exchange-correlation”  paxt  of  the  action  functional  is  for¬ 
mally  defined  as 


t  This  assumption  is  analogous  to  the  familiar  assumption  of  simulta¬ 
neous  interacting  and  non-interacting  v-representability  made  in  the 
ground-state  formalism. 
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Axc[p\  =  B,[p\  -  dtfdPrfdPr1  -  B[p]. 

Inserting  (42)  in  the  variational  equation  (33)  one  finds 


Wi  >M+^V 


*A*M.  |  =  o. 


Mr*) 


(43) 


(44) 


Since  this  equation  is  solved  only  by  the  exact  interacting  density 
one  obtains,  by  comparison  with  Eq.  (41),  the  explicit  representation 

vs(rt)  =  v(rt)+fd3r'  •  (45) 

Eqs.  (39),  (40)  and  (45)  constitute  the  time-dependent  KS  scheme. 
Compared  to  the  time-dependent  Hartree-fock  method,  this  scheme 
has  two  advantages: 

a)  The  time-dependent  effective  potential  va(r t)  is  local  (i.e.  a  mul¬ 
tiplication  operator  in  real  space). 

b)  Correlation  effects  are  included. 

Explicit  approximations  for  the  time-dependent  exchange-correlation 
potential 


(46) 

will  be  discussed  in  section  3.2  in  the  context  of  linear  response  theory. 

The  density  functionals  P[p\,  B[p\,  AXc[p\t  «tc,  are  well-defined 
only  for  “v-representable”  densities,  i.e.,  for  densities  that  come  from 
some  time-dependent  potential  satisfying  Eqs.  (4)  and  (30).  Since 
the  restriction  to  v-representable  densities  is  difficult  to  implement  in 
the  basic  variational  principle  of  section  2.2,  a  Levy-Lieb-type58,61,62 
extension  of  the  respective  functionals  to  arbitrary  (non-negative,  nor¬ 
malizable)  functions  p(rt)  appears  desirable.  Two  different  proposals 
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of  this  type  have  been  put  forward  t-j  far.63,64 

Besides  these  mathematical  generalizations,  a  number  of  exten¬ 
sions  of  the  time-dependent  density-functional  formalism  to  physically 
different  situations  have  been  developed.  Those  include  spin-polarized 
systems,65  multicomponent  systems,66  time-dependent  ensembles,67,68 
and  external  vector  potentials.64,69 

3.  FREQUENCY-DEPENDENT  LINEAR  RESPONSE 
3.1.  Selfconsistent  equations  for  the  linear  density  response 

In  this  section  we  shall  consider  the  ground-state  response  of  elec¬ 
tronic  systems.  In  that  case,  the  initial  density  Po(r)  resulting  from 
the  potential  t»o(r)  can  be  calculated  from  the  ordinary  (ground-state) 
KS  scheme: 

( -  ^lr +«0  (r)+fd3r'  +vxc\po}  (r))y?  j0)  (r)  =eipj0)(r)  (47) 

Mr)  =  £»(<,«)  Ml2-  (48) 

At  t  =  tq  a  time-dependent  perturbation  vi(rt)  is  switched  on,  i.e.  the 
total  external  potential  is  given  by 

Our  goal  is  to  calculate  the  linear  density  response  Pi(rt)  which  is 
conventionally  expressed  in  terms  of  the  full  response  function  x  48 

oi(rl)  =  Jd*  r'  J“  dt'x(Tt,  r'l')n(r'l').  (SO) 

Since  the  time-dependent  KS  equations  (39),  (40),  (45)  provide  a  for¬ 
mally  exact  way  of  calculating  the  time-dependent  density,  the  exact 
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linear  density  response  of  the  interacting  system  can  alternatively  be 
calculated  as  the  density  response  of  the  non-interacting  KS  system: 

Pi(rt)  =  f<i3r'tof00<it'xKs(rt,r,t')vl1)(^/t,).  (51) 

vi^(rt)  is  the  time-dependent  KS  potential  (45),  calculated  to  first 
order  in  the  perturbing  potential.  It  can  be  written  in  the  form 

t41}(rt)  =  vi(rt)  +  fd3r' 

+Jd3rl  Jdt1  fxe(Tt,T,t,)pi(rltl).  (52) 

fxe  denotes  the  exchange-correlation  response  kernel  which  is  formally 
defined  as  the  functional  derivative  of  the  time-dependent  exchange- 
correlation  potential  (46), 

f_r (rt  r’t')  -  (53^ 

Jxc\  rr,rtj-  gp^i^  p  =  p0,  vMJ 

evaluated  at  the  initial  ground-state  density  ?o(r). 

Eqs.  (51)  and  (52)  constitute  the  KS  equations  for  the  linear  den¬ 
sity  response.  Given  some  approximation  for  the  exchange-correlation 
kernel  fxe,  these  equations  provide  a  selfconsistent  scheme  to  calculate 
the  density  response  p\  (rt) . 

The  response  function  XKS  **  relatively  easy  to  calculate35,70  (as 
opposed  to  the  full  response  function  x)-  In  terms  of  the  stationary 
KS  orbitals  of  Eq.  (47),  the  Fourier  transform  of  XKS  (rt,  T>tl)  with  *e- 
spect  to  (t— t')  is  given  by 


XKs(ry;u)  =  Ekj(fk-fj) 


^i0)  (r)*40)  (tM0)  w *^i0)  (r') 


«  -  (ej  -  ek)  +  *6 


where  /*  and  /y  are  Fermi  occupation  factors.  The  summations  in 
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(54)  run  over  all  KS  orbitals  (including  the  continuum  states). 

For  the  purpose  of  constructing  approximations  to  the  exchange- 
correlation  kernel  fxe,  it  is  useful  to  express  this  quantity  in  terms  of 
the  full  response  function  x •  A  formally  exact  representation  of  fxc  is 
readily  obtained  by  solving  Eq.  (50)  for  vi  and  inserting  the  result  in 
(52).  Eq.  (51)  then  yields 

/*c(rt,rV)  =  x*s  (rf ,  rV) — x  “ 1  (rt , rV)  —  •  (55) 

By  virtue  of  the  Runge-Gross  theorem,  two  potentials 

u(rt)  =  t/o(r)+vi(rt)  and  v'(rt)  =  vo(r)+Vi(rt) 

always  lead  to  different  densities  p(rt)  and  p'{rt)  provided  that  vx(rt)# 
uj(rt)  +  c(t).  This  statement  does  not  automatically  imply  that  the 
linear  density  responses  Pi(rt)  and  p\ (rt)  produced  by  «i(rt)  and 
v^rt)  are  different;  the  difference  between  p(rt)  and  p'(rt)  might  show 
up  only  in  higher  perturbative  orders.  However,  since  the  right-hand 
side  of  Eq.  (19)  is  manifestly  of  first  order  in  the  perturbation,  p\{rt) 
and  Pi(rt)  must  in  fact  be  different.  Therefore  the  time-dependent 
KS  scheme  for  the  linear  density  response  (Eqs.  (51)  and  (52))  is 
well-defined  and  the  inverse  response  operators  x-1  and  x](s  *n 
(55)  must  exist71  on  the  set  of  v-representable  densities  that  corre¬ 
spond  to  the  set  of  external  potentials  specified  in  Sec.  2.1.  On 
a  larger  set  of  potentials,  the  response  functions  x  and  XKS  are 
not  necessarily  invertible.  In  particular,  the  frequency-dependent  re¬ 
sponse  operators  x(r,r/;w)  and  Xffs(r. r'jw)  cannot  always  be  in¬ 
verted  since  they  determine,  by  definition,  the  linear  density  response 
to  adiabaticallv  switched  monochromatic  perturbations, 

vj(rt)  =  vi(rw)  e~*wt  ,  -oo<f<oo, 
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which  are  not  covered  by  the  Runge-Groas  theorem  (due  to  the  condi¬ 
tion  (4)).  As  a  matter  of  fact,  it  has  been  demonstrated72,73  that  the 
response  function  x^r'ju;)  of  suitably  constructed  non-interacting 
systems  can  have  eigenvectors  with  vanishing  eigenvalues  at  isolated 
frequencies  ujq.  Consequently  the  adiabatic  density  response  vanishes 
at  these  frequencies  and  x(r,  r'\  wq)  is  not  invertible. 

We  conclude  this  section  with  a  discussion  of  the  homogeneous  elec¬ 
tron  gas,  perturbed  by  an  adiabatically  switched  potential  consisting 
of  a  single  Fourier  component  in  (gw)-space: 

t>i(rt)  =  e*(<lr-wt)Vl(qW)  e-nl*|f  —oo<t<ao. 

Again  the  Runge-Gross  theorem  is  not  applicable  in  this  case. 
However,  for  this  particular  system  the  exchange-correlation  kernel 
can  be  constructed  in  explicit  terms:  The  linear  density  response 
p1(qw)  resulting  from  the  perturbation  Vi(qw)  is  determined  by  the 
full  response  function  x^om(qw)  of  the  homogeneous  electron  gas, 

pi(qw)  =  (56) 


Comparison  of  Eq.  (58)  with  Eqs.  (51)-(52)  shows  that  the  Fourier 
transform  /£°m(r  -  r\t  -  t1)  of  f^°m(qu)  is  indeed  the  exchange- 
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correlation  kernel  of  the  homogeneous  electron  gas.  /«m(?w)  is 
closely  related  to  the  so-called  local  field  correction  G(qu)  of  the  ho¬ 
mogeneous  electron  gas: 

(59) 

Using  well-known  features  of  the  electron  gas  response  functions,  f£°' m 
can  be  shown  to  have  the  following  exact  properties: 

1.  As  a  consequence  of  the  compressibility  sum  rule  one  finds74 

lim /£com(g,w  =  0)  =  j^(pe*c(p))  =  fo{p),  (60) 

g— *  0  dp 

where  exc(p)  denotes  the  exchange-correlation  energy  per  particle  of 
the  homogeneous  electron  gas. 

2.  The  third-frequency-moment  sum  rule  leads  to75 

Km w  =  oc) 

<7—0 

=  s  /-«•  <6l> 

3.  According  to  the  best  estimates78,77  of  exc(p),  the  following  relation 
holds  for  all  densities 

fo(p)  <  foo(p)  <  0.  (62) 

4.  In  the  static  limit  (w  =  0),  the  short-wavelength  behavior  is  given 
by78 
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?Um  /*«”(,,„  =0)  =  (1  -  j(0)),  (63) 

where  y(r)  denotes  the  pair  correlation  function. 

5.  For  frequencies  w^O,  the  short-wavelength  behavior  is79 

?Km  /£?m(4.^0)  =  -§  lf(l- ,(<>)).  (64) 

6.  fx?m{qu>)  is  a  complex- valued  function  satisfying  the  symmetry  re¬ 
lations 


Ref^r(Q^)  =  -")  (65) 

Im  /*?"»(*,  ")  =  -/m  &om(q,  -w).  (66) 

7.  f^°m(qu)  is  an  analytic  function  of  w  in  the  upper  half  of  the  com¬ 
plex  w-plane  and  approaches  a  real  function  /<»(?)  for  u;— +oo.80  There¬ 
fore,  the  function  /*°m(9<*>)  — /oo(?)  satisfies  standard  Kramers-Kronig 
relations: 

/*r(?w)  -  /«(,)  =  p/r^  ^  («») 

imftrfa)  =  -p/?„  ¥  Rt  Mq)  ■  (68> 

8.  The  imaginary  part  of  /££"*  exhibits  the  high-frequency  behavior 

JimlmfirM  =  3/2  (6») 

for  any  q  <  oo.81  A  second-order  perturbation  expansion81,82  of  the 
irreducible  polarization  propagator  leads  to  the  high-density  limit 
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c  =  23W15.  (70) 

9.  In  the  same  limit,  the  real  part  of  /£com  behaves  like83 

lim  RefiT(q, w)  =  /oo(g)  +  -fa.  (71) 

Since  c  >  0,  the  infinite-frequency  value  i»  approached  from  above. 
This  implies,  in  view  of  the  relation  (62),  that  Re  f£?m(q  =  0,  w)  can¬ 
not  grow  monotonically  from  /o  to  /oo- 

The  above  features  of  fx°m  are  valid  for  a  three-dimensional  elec¬ 
tron  gas.  Analogous  results  have  been  obtained  for  the  two-dimensional 


3.2.  Approximations  for  the  exchange-correlation  kernel 

The  time-dependent  density-functional  formalism  described  in  sec¬ 
tion  2  can  be  shown51  to  reduce  to  the  traditional  HKS  formalism  for 
the  class  of  static  initial  value  problems  in  which  the  potential  is  time- 
independent,  v(rt)  =  v(r),  and  where  the  initial  state  ‘s  the  ground 
state  of  v(r).  In  particular,  the  functional  AXc[p\  reduces  to  the  usual 
exchange-correlation  energy  if  ground-state  densities  are  inserted  for 
p(rt): 

AZe[P\  =  (ti-to)'£xc[0]  for  all  ground -state  densities p(r).  (72) 

This  has  an  immediate  consequence  for  possible  approximations  to 
the  functional  ^lxc[p]*  1^  an  approximate  expression  for  ^xc  \p\  both 
satisfies  (72)  and  is  local  in  time,  it  must  necessarily  have  the  form 


(73) 
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Of  particular  interest  is  the  adiabatic  local  density  approximation 


A*LDA[p\  =  f^dtfd3r  [p-€xe{p)\ |p=p(rt)  (74) 

which  is  local  in  time  and  space.  By  Eqs.  (46)  and  (53),  the  exchange- 
correlation  kernel  corresponding  to  (74)  takes  the  form 


f£LDA(Ttyt')  = 


(75) 


Obviously,  the  static  exchange-correlation  energy  can  be  expected  to 
be  a  good  approximation  only  for  very  slow  time-dependent  processes, 
i.e.,  one  implicitly  invokes  an  adiabatic  approximation.  This  is  also 
reflected  in  the  fact  that  (75)  is  identical  with  the  zero-freouencv  limit 
(60)  of  the  homogeneous  exchange-correlation  kernel: 


tALDA 
J  xe 


(r  y-.w)  =  5(r-rO/rCom(9  =  0,w  =  0;p) 

0= 


P=Po(r) 


(76) 


Going  beyond  the  adiabatic  local  density  approximation,  Gross  and 
Kohn83  have  proposed  an  approximation  for  the  exchange-correlation 
kernel  which  explicitly  incorporates  the  frequency  dependence.  With 
respect  to  spatial  coordinates,  a  local  approximation  as  in  (76)  is  em¬ 
ployed: 

=  «(r-r')/ir(«=<W)|  (77) 

>P=Po(r) 

The  frequency-dependence  of  I m  /£com  (g  =  0,o»)  is  approximated  by 
a  Pad^-like  interpolation  between  the  high-  and  the  low-frequency  lim¬ 
its: 


=  O.a,) 


o(p)‘U> 


(78) 
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where 

«(p)  =  -<W4s,s(/«W-/oW)s^ 

=  b/«)4^(/oo((.)-/0(0))4''5 

/Oi/oo»  “d  c  are  given  by  Eqs.  (60),  (61),  and  (70),  respectively. 
For  the  exchange-correlation  energy  ezc(p)  (which  enters  /o  and  /oo)> 
the  parametrization  of  Vosko,  Wilk  and  Nusair77  is  employed.  Using 
the  Kramers-Kronig  relation  (67),  the  real  part  can  be  expressed  as 

=  /oo  +  2£d/V2)  -  4^n(l^,-^) 

-  •>2  =  i  +  ‘“2-  P») 

E  and  II  are  complete  elliptic  integrals  of  the  second  and  third  kind  in 
the  standard  notation  of  Byrd  and  Friedman.86 

Figs.  1  and  2  show  the  real  and  imaginary  part  of  fx°m  as  cal¬ 
culated  from  (78)  and  (79).  The  functions  aire  plotted  for  the  two 
density  values  corresponding  to  ra  =  2  and  rs  =  4.  For  the  lower  den¬ 
sity  value  (r,  =  4),  a  considerable  frequency-dependence  is  found.  The 
dependence  on  ui  becomes  less  pronounced  for  higher  densities.  In  the 
extreme  high-density  limit,  the  difference  between  /o  and  /<»  tends  to 
zero.  One  finds  the  exact  result 

foo-fo  ~  r\ 


for  r,— »0. 


(80) 


Fig. 


fuj(a.u.) 


1.  Real  part  of  the  parametrization  for  /jCom(g  =  0,w) 
(from  Ref.  75). 


Fig.  2.  Imaginary  part  of  the  parametrization  for 
fxcm (f  =  0, u/)  (from  Ref.  75) . 
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At  the  s&me  time,  the  depth  of  the  minimum  of  Imf^m  decreases, 
again  proportional  to  r2. 

We  finally  mention  that  an  extension  of  the  parametrization  (78)  to 
nonvanishing  q  was  given  by  Dabrowski.87  A  similar  interpolation  for 
the  exchange-correlation  kernel  of  the  2-dimensional  electron  gas  has 
been  derived  by  Holas  and  Singwi.81 

3.3.  Applications 

The  time-dependent  KS  scheme  defined  by  Eqs.  (51)  and  (52)  has 
been  successfully  applied  to  the  photo-response  of  atoms35,88,89  and 
molecules,90,91  metallic92-98  and  semiconductor99  surfaces  and,  most 
recently,  of  bulk  semiconductors.100 

In  the  present  context,  we  restrict  ourselves  to  sufficiently  low  ra¬ 
diation  frequencies  such  that  the  wavelength  is  much  larger  than  the 
atomic  radius  and  that  Compton  scattering  is  negligible.  In  the  case 
of  atoms,  this  condition  is  satisfied  for  photon  energies  below  «3  ktV . 
Under  these  circumstances,  the  electric  field  across  the  atom  can  be  as¬ 
sumed  to  be  constant. 

The  potential  vi(rt)  associated  with  a  monochromatic  spatially 
constant  electric  field  in  z-direction  is  given  by 

vj(r<)  =  Ezcoeut.  (81) 

The  induced  density  change,  6p(rt),  is  most  conveniently  character¬ 
ized  by  a  series  of  multipole  moments.  The  induced  dipole  moment, 

P(0  =  -/ zSp(rt)d3r,  (82) 

can  be  expanded  as101 


p(t )  =  a(w)-Ecoswt+J/J(0)  :  EE+J/?(2w)  :  EE  cos  2  wt 
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+  |'r(w)  EEEcosu;t+|'7(3w)  :  EEEcoe3wt,  (83) 

where  the  notation  is  meant  to  convey  the  tenaorial  character  of  the 
quantities  a,/3,7.  The  first  coefficient,  a,  is  termed  the  (dipole)  po¬ 
larizability,  0  and  "  denote  the  second-and  third-order  ( dipole )  hyper¬ 
polarizabilities.  For  ordinary  laboratory  fields,  the  electric  dipole  mo¬ 
ment  is  adequately  described  by  the  polarizability  a.  In  the  presence 
of  intense  laser  fields,  hyperpolarisabilities  up  to  7  have  to  be  taken 
into  account.  Higher  orders  are  negligible  in  present-day  experiments. 
For  spherically  symmetric  states  0  is  zero. 

If  one  considers  only  the  linear  term  in  Eq.  (83)),  the  Fourier  com¬ 
ponents  Pi(rw)  of  the  first-order  density  shift  Pi(rt)  can  be  treated 
independently  and  the  frequency-dependent  polarizability  is  given  by 

=  -%fzPl(™)d3r.  (84) 

By  virtue  of  the  Golden  Rule,  the  photoabsorption  cross  section  is 
then  proportional  to  the  imaginary  part  of  a(w): 

c(u)  =  Ima(v).  (85) 

Zangwill  and  Soven35  have  used  the  time-dependent  KS  scheme  (51)- 
(52)  to  calculate  the  linear  density  response  of  rare-gas  atoms  to  a 
perturbation  of  the  form  (81).  For  the  exchange-correlation  kernel 
fxc  they  employed  the  adiabatic  local  density  approximation  (76). 
The  resulting  photoabsorption  cross  section  (as  calculated  from  Eqs. 
(84)  and  (85))  turns  out  to  be  in  rather  good  agreement  with  exper¬ 
iment.  As  an  example,  we  show  in  Fig.  3  the  photoabsorption  cross 
section  of  the  Xe  atom  in  the  vicinity  of  the  4d  threshold.  A  few 
general  remarks  are  in  order  at  this  point: 


6  7  8  9 
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Fig.  3.  Total  photoabsorption  cross  section  of  the  Xe  atom  versus 
photon  energy  in  the  vicinity  of  the  4d  threshold  (from  Ref.  35).  Solid 
line:  selfconsistent  time-dependent  KS  calculation;  dash-dot  line:  self- 
consistent  time-dependent  Hartree  calculation  (/zc  =  0);  dashed  line: 
independent  particle  result  (Hartree  and  exchange-correlation  kernels 
neglected);  crosses:  experimental  data  from  Ref.  102. 

i)The  calculations  of  Zangwill  and  Soven  do  not  show  the  char¬ 
acteristic  autoionization  resonances  found  in  high-precision  measure¬ 
ments.  This  is  due  to  the  local  density  approximation  employed  for 
the  ground-state  exchange-correlation  potential  in  Eq.  (47).  Since 
self-Coulomb  and  self-exchange  energies  do  not  properly  cancel  in  the 
local  density  approximation,  the  total  selfconsistent  potential  in  Eq. 
(47)  falls  off  exponentially  rather  than  like  -1/r.  As  a  consequence, 
the  effective  potential  does  not  support  the  high-lying  Rydberg  states 
essential  for  the  description  of  autoionization  resonances.  In  princi¬ 
ple,  the  time-dependent  KS  scheme  (51)-(52)  should  be  capable  of 
describing  such  effects  if  a  selfinteraction-corrected  approximation  for 
the  exchange-correlation  potential  is  used.  A  different  approach  to  the 
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selfinteraction  problem  was  proposed  by  Zangwill  and  Soven  in  a  sec¬ 
ond  paper,89  where  they  introduced  selfenergy  corrections  directly  in 
the  response  function  XKS- 

ii)  Within  the  original  scheme  of  Zangwill  and  Soven,  the  onset  of 
continuous  absorption  for  each  atomic  subshell  is  given  by  the  corre¬ 
sponding  KS  single-particle  energy  eigenvalue.  This  follows  directly 
from  the  fact  that  the  response  function  XKS  has  poles  at  the  KS 
single-particle  excitation  energies  (cf.  Eq.  (54)).  Since  the  KS  ex¬ 
citation  energies  do  not  coincide  with  the  exact  atomic  excitation 
energies,  the  calculated  absorption  edges  are  typically  several  volts 
below  the  observed  thresholds.  We  emphasize  that  this  result  is  due 
to  the  approximation  (76)  used  for  /xc(r,  r'lu/).  The  exact  function 
fxc  (r-r/;u),  formally  defined  as  the  Fourier  transform  of  (55),  has  a 
rather  intricate  structure  which  both  removes  the  incorrect  poles  of 
XKS  gives  rise  to  poles  at  the  correct  quasi-particle  excitation  en¬ 
ergies. 

iii)  It  is  of  interest  to  estimate  the  overall  importance  of  the  exchange- 
correlation  kernel  /zc.  The  dash-dot  curve  in  Fig.  3  shows  the  result 
obtained  with  fxe  set  equal  to  zero.  Deviations  from  the  full  calcula¬ 
tion  (solid  curve)  are  about  10%.  Clearly  the  full  calculation  agrees 
much  better  with  the  experimental  values. 

iv)  Zangwill  and  Soven  have  employed  the  adiabatic  local  density 
approximation  (76),  in  which  /*e(r,r,;w)  is  replaced  by  /*c(r,  r7;  w  = 
0).  On  the  basis  of  the  parametrization  (78)-(79)  of  Gross  and  Kohn,83 
the  importance  of  the  frequency  dependence  of  the  exchange-correlation 
kernel  has  been  estimated:  In  the  frequency  range  relevant  for  atomic 
photoabsorption,  the  ^-dependent  function  fXc  is  expected75  to  differ 
by  at  most  2-6%  from  its  zero-frequency  value  employed  in  the  calcula¬ 
tion  of  Zangwill  and  Soven.  Since  the  complete  neglect  of  fxc  leads  to 
deviations  of  only  10%,  the  effect  of  the  frequency-dependence  of  fxc 
can  be  estimated  to  be  smaller  than  1%  in  the  photoabsorption  cross 
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section  of  atoms. 

The  scheme  of  Zangwill  and  Soven  has  also  been  applied  to  calcu¬ 
lations  of  the  photoabsorption  of  small  molecules.90,91  In  particular,  a 
previously  unexplained  double-peak  structure  near  15  eV  in  the  pho¬ 
toemission  cross  section  of  acetylene  could  be  described  successfully. 

The  adiabatic  exchange-correlation  functional  (74)  has  also  been 
employed  in  calculations  of  the  third-order  hyperpolarizabilities  de¬ 
fined  in  Eq.  (83).  The  results  obtained  by  Zangwill103  for  rare-gas 
atoms  are  in  good  agreement  with  experiment.  According  to  Senatore 
and  Subbaswamy,104  however,  a  term  was  missing  in  Zangwill’s  calcu¬ 
lation.  When  they  include  this  term,  the  agreement  with  experiment 
is  reduced  considerably:  The  calculated  third-order  hyperpolarizabil¬ 
ities  of  rare-gas  atoms  are  too  large  typically  by  a  factor  of  2. 

The  response  problem  of  surface  electrons  is  to  find  the  currents  and 
charge  densities  induced  by  perturbing  electromagnetic  fields,  and  to 
calculate  the  resulting  modifications  of  these  fields. 

The  Hamiltonian  governing  the  electronic  dynamics  is 

^ =E,  j  (Pi  ~  ?  A  x  (^ ,  t) ) 2 + Ei  ( vo  (ri ) + vi  (^ ,  t) ) 

+  2  |r,  -rj|’  <86) 

where  t>o(r)  i®  the  potential  due  to  the  ion  cores  and  vi(r,t),  Ai(r,t) 
are  the  perturbing  potentials  in  the  transverse  gauge  (V-Ai  =  0).  The 
primary  objective  is  to  find  the  first-order  density  and  current  changes, 
Pi(r,t)  andji(r,t),  caused  by  Vi  and  Aj. 

In  most  theoretical  work,  as  well  as  in  this  section,  the  ionic  po¬ 
tential,  vq(r),  is  replaced  by  the  potential  due  to  a  uniform  positive 
charge  background  in  a  half-space,  say  z  >  0.  This  is  the  so-called 
jellium  model  for  metal  surfaces.  In  this  model  there  are  two  intrinsic 
microscopic  length  scales,  the  inverse  Fermi  wave-number,  kjl,  and 
the  Thomas-Fermi  screening  length  («  surface  thickness),  In  the 
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present  context,  both  lengths  are  typically  of  the  order  a«10  8  cm. 

In  the  most  important  applications  the  perturbing  potentials  vj 
and  A}  vary  on  a  length  scale,  l,  which  satisfies  l  a.  Examples 
are  the  scalar  potential  vi(r)  due  to  an  external  charge  at  a  distance 
z'Sxi,  or  the  vector  potential,  associated  with  a  light  wave  of 

wavelength  A »a.  The  corresponding  ft*  and  j\  vary  on  the  scale  of  l 
in  the  x-y  plane  but,  because  of  the  abrupt  drop  of  the  unperturbed 
density  at  the  surface  (on  the  scale  of  a),  they  vary  on  the  short 
scale  a  in  the  z-direction.  Formal  arguments  due  to  Feibelman105 
have  shown  that,  to  leading  order  in  a/ty  the  effect  of  the  surf  _e 
on  the  electromagnetic  fields  far  from  the  surface  (|z|  >  a)  is  entirely 
characterized  by  two  complex  frequency-dependent  lengths,  d|j(w)  and 
d±{uj). 

Actually,  for  the  jellium  model  d||(w)=0.*  To  define  dj.(w)  we 
Fourier  analyze  all  physical  quantities  parallel  to  the  surface,  in  the  x- 
y  plane.  For  example,  a  Fourier  component  of  the  induced  charge  den¬ 
sity  becomes 

Pi(r,w)  =  (87) 

where  qj|  =  {qx,qy, 0)  (and  |qy|-a<l).  Then  djjw)  is  given  by 

=  fzpi{z,v)dz/fpi{z,w)dz,  (88) 

i.e.  it  is  the  (complex)  center  of  mass  of  the  induced  surface  charge. 
d±(u)  is  the  generalization  of  the  static  image  plane  introduced  by 
Lang  and  Kohn.10® 

This  di(w)  is  then  the  object  of  quantitative  calculations.  They 

*  This  result  has  been  obtained  in  the  random  phase  approximation  in 
Ref.  105.  It  is  easily  established  as  a  rigorous  many-body  result  for 
the  jellium  model  (W.  Kohn,  unpublished). 
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require  the  density  response,  pj(z,w),  to  a  uniform  exte..ial  electric 
field  perpendicular  to  the  surface.  The  calculation  was  first  carried 
out  in  the  random  phase  approximation  (RPA),  equivalent  to  time- 
dependent  Hartree  theory,  in  which  effects  of  exchange  and  correlation 
are  neglected  in  the  self  consistent  perturbing  potential  (52).  These 
calculations  led  to  very  interesting  results  not  present  in  classical 
Maxwell  theory,  the  surface  photo  effect,  excitations  of  bulk  and  sur¬ 
face  plasmons,  and  hints  of  an  additional  surface-localized  plasmon. 
They  help  to  explain  experimental  data  on  the  photo-yield  spectrum 
of  At 107 

Some  recent  calculations92-98  have  employed  the  time-dependent 
response  theory  of  Gross  and  Kohn,  including  exchange-correlation  ef¬ 
fects  either  within  the  adiabatic  approximation  (75)  or  the  frequency- 
dependent  parametrization  (78)-(79).  They  have  confirmed  the  qual¬ 
itative  conclusions  of  the  RPA  calculations,  with  quantitative  modifi¬ 
cations  generally  of  the  order  of  10%. 

Finally  we  mention  some  recent  formal  work  on  the  response  prop¬ 
erties  of  real  crystal  surfaces  (not  the  jellium  model).108  Here  again 
the  response  properties  of  the  surface  are  characterized  by  complex 
d- parameters  having  the  dimension  of  a  length.  Quantitative  calcula¬ 
tions  have  been  reported108  to  be  in  progress. 
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